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1 Introduction 



The purpose of this paper is to make a detailed analysis of the one-point functions in a 
particular integrable boundary quantum field theory, the boundary scaling Lee- Yang model. 
Such a study is of intrinsic interest, allowing one to check the consistency of various different 
approaches to the computation of these quantities. It also serves to illuminate some issues 
that have arisen in recent investigations of the spectra and boundary entropies ((7~functions) 
of integrable models with boundaries. 

After a brief review of the pertinent features of the boundary scaling Lee- Yang model in 
section 2, sections 3 and 4 of the paper are devoted to the one-point functions of bulk fields 
in the presence of a boundary. Since these generally have a dependence on the distance from 
the boundary, their exact determination is a non-trivial problem, of similar difficulty to the 
calculation of two-point functions in the absence of boundaries. For this reason, apart from 
some simple cases associated with the Ising model |13], no exact formulae for these boundary 
one-point functions are known. Instead, one has to resort to approximate methods, two of 
which we investigate in some detail. The first is the truncated conformal space approach, or 
TCSA. This has previously been used to study integrable models both without [@|2^ and 
with 1^,^ boundaries, but this particular application of the technique is new. The second 
method is based on a form factor (FF) expansion, and makes use of the expansion of the 
boundary states in the basis of infrared multi-particle states given in |11|. Previous work on 
this topic includes refs. [13|, 14, 17, 1^]. The main novelties arising in the discussion of the 
scaling Lee- Yang model presented below are that this model has a non-trivial bulk S-matrix, 
and that the boundary one-point functions receive contributions from both even- and odd- 
particle-number components of the boundary state. This allows the relative normalisations of 
the two sectors to be checked, and we find evidence that the prescription given in |11| is out by 
a factor of 2. Section 3 introduces and discusses the TCSA and FF methods, culminating in 
a numerical comparison between the two which, modulo the small correction to the boundary 
state just mentioned, shows an excellent agreement. 

Both the TCSA and the FF methods can be pushed a little further, allowing expectation 
values in states other than the ground state to be accessed. These generalisations are explained 
in section 4. 

The expectation values of boundary fields are also of interest, and for these there is more 
hope to find exact results, at least for certain geometries. This is the subject of section |^. 
First, in section 5.1, we describe the so-called 'R-channel' approach, which allows us to obtain 
the expectation values of boundary fields at the edge of a strip of finite width, both in the 
ground state and, in principle, in any excited state. By relating the expectation values to 
the derivative of the energy with respect to the boundary field, we find that they can be 
expressed in terms of the solutions to TBA equations, and compare these results with data 
from the TCSA (obtained by adapting the method of |12|), for various different combinations 
of boundary conditions and for varying strip width. Then in section 5.2 we switch to the 
'L-channel'. Making use of results from we are able to relate the expectation value of a 
boundary field placed at one end of a semi-infinite cylinder to the so-called 'Y-function' of the 
Lee- Yang model. An interesting relationship between the spectra of certain different models 
placed on the same strip emerges as a by-product of this discussion. 

Finally, in section ^, we make use of the results in the preceding sections to examine 
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the RG flow from tlie ^{h) boundary to the 1 boundary, first discussed in and to treat 
some previously obscure features of the model on a strip of finite width, with simultaneous 
perturbations on both boundaries. Section |^ contains our conclusions, and indicates some 
directions for future work. 

As in previous papers ||5|,|^|8|, we will be using the boundary scaling Lee- Yang model as 
our example. We have included a brief review of the features of this model in the next section, 
but the earlier works should be consulted for some more detailed explanations. 

2 The boundary Lee- Yang model 

2.1 The conformal field theory description of the critical Lee- Yang model 

The Lee- Yang model is the simplest non-unitary conformal field theory, A^2,5, and has central 
charge —22/5 and effective central charge 2/5. There are only two representations of the 
Virasoro algebra of interest, of weight and —1/5, and consequently only two bulk primary 
fields, the identity 1 of weight 0, and ip of weight x^ = — 2/5; equally there are only two 
conformally- invariant boundary conditions, which we denote by 1 and ^. 

There are three non-trivial boundary fields interpolating the various boundary con- 

ditions a and /?, of weight = 1/5. Two of these {1^,1^^) interpolate the two different 
conformal boundary conditions^, and one {(f>) lives on the ^ boundary: 

(2.1) 

The OPEs of interest are the bulk OPE, 

ip{z,z) ip{w,w) = C^^^\z-w\^^^ + C^^"-^ \z - wl"^^^ ip{w,w) + ... , 

the boundary OPEs, 

Hz)^{w) = C4,/\z-w\^/'' + C^/\z-w\^/''(l)iw) + ... , 

^l^{z)4>{w) = C^/\z-w\^/^^l;{w) + ... , 
^{z)i^Hw) = C^^f\z-w\'/'i,Hw) + ... , 
^(z)^t(^) = C^^,^\z-w\^/' + ... , 

^Hz)i^{w) = C^,/\z-wf^' + C^t/\z - w\^/'' cPiw) + ... , 
and the two bulk-boundary OPEs 

= (^)i?^|2(z-u;)|2/5 + ... , 
(p(z)l^ = (*W|2(z-«;)|2/5 + (*'i?$|2(z-^i;)|i/5</,(u;) + ... . 



*In M we did not distinguish these fields and denoted them both by ip 
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A suitable choice for the structure constants is 
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There are three possible choices for pairs of conformal boundary conditions on a strip: (1,1), 
($, 1) and (<I>, <!>). We shall take the strip to be of width R with coordinates < x < across 
the strip and y running along the strip, and normalise all our correlation functions so that 
the expectation value of the identity operator on a strip is always one. 

The strip correlation functions between states ( a \ and | (3 ) can be found by mapping the 
strip to the unit disc and inserting the appropriate fields il^a and ijjp. Since the ground state 
on the strips with boundary conditions (1, 1), (1,*!*) and ($,<!?) correspond to the fields 1, 
and (f) respectively, one needs to include the appropriate field insertions to find the ground 
state expectation values on these strips. These insertions lead directly to the particular chiral 
blocks and structure constants in (|2.4|) and (^). 

The one-point functions of the field (/?(x) on such a strip are best expressed in terms of 
the four strip chiral block functions fi{0), 
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In terms of these functions, we have 
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It is only for the latter two pairs of boundary conditions that the boundary field (p{y) exists 
on the boundary x = 0. For these cases the one-point functions are simply 
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Of the five expectation values ( |2.4D , (|2.5| ) , only the first has a finite limit as the strip width 
tends to oo, the others all diverging. None of them depend on y, and so when no confusion 
can arise this variable will often be omitted, even inside the vacuum expectation values. 



2.2 The scaling Lee- Yang model 

The scaling Lee- Yang (SLY) model can be described as a perturbation of the critical Lee- Yang 
model by the term 

A / Lp{w,w) d^w . (2.6) 



This leads to a massive scattering theory, comprising a single particle with two-particle S 
matrix 

sinh(| + ^) 



S{e) = -{1){2) , (x) 



sinh 



'8 _ iTTx\ ■ 
.2 6 / 

The mass M of the particle is related to the bulk perturbation parameter A by |p4, 26| 



M = ACA5/12 



(r(3/5)r(4/5))^/^2 



2.642944... . 



(2.7) 
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55/i6r(2/3)r(5/6) 

We will also need the form factors of the bulk model. These are the matrix elements of 
the elementary field (p{x) in the asymptotic n-particle states which can be formally written 
in terms of the ZF operators A{9) as \ 9i, . . . ,6^) = ^(^i) ... A{9n)\0). The form-factor 
Fn{9i ...9n) is then given by 

(O|(^(x,O)|0i,...,0„) =exp(-MxE^cosh0,)^«(^i,---,^n) • (2.9) 

The form factors of the SLY model were first computed in |2l|; we, however, adopt the 
conventions of [25|, modulo the fact that for us ip{x) is a real field. The function F„ can be 
parametrised as [25|: 



1 ^n) — HnQni^i, 



n 

,Xn)Y[ 

i<j 



Xi ~\~ X 



(2.10) 
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where Xi = exp(0j), i = 1 
factor bootstrap 



The terms in ( 2.10 ) can be determined through the form 



with the result 



v{i7r — 9) v{—iTT + 



cosh 6* + 1/2 

where we take the function in a form suitable for numerical evaluation |2^ 

N 



(2.11) 



v{9) 



n 

n=l 



+ n + l/2)(^ +n- l/6)(^ +n- 1/3) 
(2^+^-V2)(o!,+n+l/6)(^+n+l/3) 



X exp 2 



dt 



sinh(t/2) sinh(t/3) sinh(t/6) 
tsinh2(t) 



(AT + 1 _ iVe-2t) ^-2Nt+i9t/n 



viO) 



1.111544045... 



(2.12) 
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(with N arbitrary) and 



_i3^ 
V2v{0) 



(2.13) 



This overall normalisation of the form factors is taken from the results of |10|, where (in our 
conventions) the expectation value ( ) in the bulk is 
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Using the relation between M and A in (|2.8| ), this boils down to 

-3^r(i)f M-f 



(2vr)V54r(i)r(|) 



(-1.239394325... )M~2/^ . 



(2.15) 



The functions Qn{xi, ■ ■ ■ ,x„) in ( 2.101) are symmetric polynomials of degree n(n — l)/2 and 
partial degree n — 1. These polynomials have been determined via the form factor bootstrap 
approach for arbitrary particle numbers n. They can be nicely written in the form of a 



determinant of a matrix in symmetric polynomials |21,25], for a related formulation see |20| 
For our purposes it will be sufficient to list the first few: 
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(2.16) 



where the elementary symmetric polynomials in n variables e)- are defined by 

n n 



i=l 



r=0 



The integrable boundary conditions for the model were discussed in detail in [^. The 
allowed boundary conditions are the 1 conformal b.c, and the perturbation ^{h) of the 
conformal <I> boundary by the integral along the boundary 

h j (p{x) dx . (2.17) 

The boundary reflection factors corresponding to these two boundary conditions are 

R^^h){0) = Rb{0) , R^{0) = Ro{e), (2.18) 

where 



MO) = a) (1) ilV' {3(9 + ^^'-^)S{e - ^n!^)y' . 
The relation between b and h was conjectured in Q to be 

Hb) = -\h^,it\ sin(^(6 + l/2)/5) . 



(2.19) 
(2.20) 
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Sometimes we shall find it useful to consider instead the dimensionless quantity 



^crit — hcvitM , 

a constant which was found in [^j: 

/.crit = -vr^/' 2^/5 5V4_j_|__ f _A|Z j = -0.68528998399118 .... (2.21) 

Combining (^]^) and ( |2.21] ), we obtain the more convenient formula 

53/8 sin(^(fe+l/2)/5) 

^-~¥J^ — ^ ■ ^^-^^^ 

Finally, we will need the boundary-particle couplings for the various boundary condi- 
tions. In |]ll[ these were defined in two different ways, either via the residue at 9 = iiT/2 of 
the reflection factor R%^{9) for a particle of type a on the boundary a : 

K{0)-^^^^ (2.23) 
2 a — ITT 1 1 

or, in models where bulk fusings occur, via the residues of certain other poles, divided by 
the corresponding bulk couplings. (If the boundary scattering is non-diagonal, the formulae 
become a little more complicated - see [^.) For the Lee- Yang model, the consistency of the 
two definitions was shown to follow from the bootstrap equations and crossing in Since 
for this case there is only one particle type, the particle index a can be dropped, and we have 

• For the 1 boundary. 



g^ = -i2\J 2\/ 'i - 3 . (2.24) 

For the ^{h[h)) boundary, 

tan((6-F2)7r/12) , ^ , 



Notice that while the reflection factors for the 1 and $(/i(0)) are identical, the corresponding 
boundary-particle couplings differ by a sign. This can be traced to the fact that when 6 = 0, 
the residue of at zvr/G for the <I>(/i(0)) boundary receives additional contributions from 
intermediate states containing boundary bound states. The net effect is to negate the cou- 
pling; indeed, this is the only option given that the residue at z7r/2 must remain unchanged. 
In terms of the boundary states to be discussed in the next section, this means that the only 
difference between the (infinite-volume) boundary states for these two boundary conditions 
is in the sign of the contributions from states of odd particle number. 
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3 The one-point functions of the bulk field 



3.1 {^(x) ) from form factors 

In this section we consider the one-point function ( ip{x) ) in the presence of a boundary using 
the form factor approach. We will restrict consideration to a theory with a single scalar 
particle, as in the Lee- Yang model. Schematically, the idea of the form-factor approach is to 
evaluate the one-point function on the upper half plane as 

oo 

{ifix)) = {0\^{x)\B^) = J2 (019^(^)1^) {n\B^), 

n=Q 

where | ) is a boundary state corresponding to the boundary condition a and the sum 
over asymptotic states has been split into the contributions from n = 0, 1, 2, ... particles 
(here \n){n\ represents the projection onto asymptotic states with n particles). We have 
also taken the boundary state to be normalised such that ( {Bq, ) = 1. 

This boundary state can be expanded in terms of multi-particle states on the infinite 
line, using the Zamolodchikov-Faddeev (ZF) operators A{6) which create single particles of 
rapidity 9. When the state contains no zero-rapidity particles, it can be written as |11] 



Bn 



exp 



d9_ 

4tt 



Ka{e) A{-e) A{9) 



(3.1) 



where Ka[9) is related to the reflection amplitude Ra{9) for the a boundary condition by 



(3.2) 



In general there may also be contributions to the boundary state involving zero-momentum 
particles, which can be associated with couplings of single bulk particles to the boundary. Up 



to the three-particle contribution, the appropriate boundary state was given in |11] as 

r-oo 
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1 + g^A{0) + 
9a A{0) 



Ait 



K^{9) A{-9) A{9) 



CO sin 

^K^{9)Ai-9)A{9) + 



and it is natural to suppose that the full expression is 



(3.3) 



exp 



/oo JQ 
—K^[9) A{-9) A{9) 



o; 



(3.4) 



In |^T|, the factor 'g^ was identified with the boundary-particle coupling Qa, as defined at 
the end of the last section. However, our numerical results (and also an examination of the 
reflection factor and boundary state for the Ising model with free boundary conditions given 
in ||ll|) cast doubt on this suggestion. As will be explained in section |3^ below, we found 
that we had rather to set 

ga = Qal'i (3.5) 
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in order to obtain a successful match with TCSA data. 

The form factors for the Lee- Yang model were given in the previous section, and substi- 
tuting these into (|3.3|) , we find the leading large-x behaviour of ( tp{x) ) is 



_ 33/4 d0 f(0)f(-e)f(-2e) (1+2 cosh 

-Ma;(l+2cosh(6»)) , 

2V2vi0r ' 8cosh0cosh(0/2)2 



(3.6) 



This will be compared with TCSA data in section 3.2 below. 



3.2 Estimating ( ip{x) ) using the TCSA 

The first step in the calculation of ( (p{x) ) using the TCSA is the numerical evaluation of the 
ground state | ) of the perturbed Hamiltonian with dimensionless strip width r = MR. 

H{R,X,hi,hr) 

TT 

where O represents the operator O on the upper half plane restricted to the conformal space 
truncated to level A^. The parameter A determines the bulk mass and we have allowed the 
possibility of boundary fields on the left and right edges of the strip, with strengths hi and 
hr- The second step is then to estimate the expectation value of the (dimensionless) operator 
in terms of the matrix elements of the operator <^ as 

( V(x) ) ~ (3 (0|^(exp^Wi^))|0) _ ^3.7) 

The state | ) can be expanded in Virasoro primary and descendent states; by repeatedly 
commuting Virasoro algebra elements through (p, the general matrix elements of if can be 
expressed in terms of the matrix elements between Virasoro primary states, and their deriva- 
tives. In the rest of this section we discuss the TCSA method in more detail; a comparison 
of the TCSA and FF results is given in section All the TCSA results in this paper were 
calculated on a workstation in Mathematica with truncation levels up to 18 and on spaces 
with up to 161 states. 
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3.2.1 The strip with (1,1) boundary conditions 

The simplest case is the strip with boundary conditions 1 on both sides. The unperturbed 
conformal field theory expectation value is given in (^^). If we denote the scaled position of 
the field by ^ and the normalised strip width by r where 

i = xM , r = MR , so that < ^ < r , 

then the TCSA estimate of the expectation value in the bulk perturbed model truncated to 
level N takes the form 

g(^-)(0-M2/^(^(x) ^ =(^sin^)'^'x;/^^^-^f^) • (3-8) 

^ ' N,r V vr r / ^-^ \ r J 

^ ' n=0 ^ ^ 

The coefficients f^'^'^ are determined by the expansion of the ground state | ) and have to 
be calculated numerically. Since the state | ) lies in the /i = representation, the matrix 
elements of (p are given in terms of the chiral block /a; furthermore, since this representation 
has a null state at level 1, one can eliminate all states containing -L_i, and so one does not 
need any terms with derivatives of this chiral block in (|3.8| ) . To compare with the form- factor 
calculation we will need to take the simultaneous limits — > oo, r ^ oo while keeping ^ 
fixed. We shall often drop the labels (iV, r) and write simply Q{f^ for the TCSA estimates. 

To show the typical behaviour of these quantities, we give the values of /A for r = 8 
and 4 < A^ < 12 in table |, and for A^ = 12 and 2 < r < 12 in table |. In figure [l| we 
plot ^^^^''"^(c^) for fixed truncation level A^ = 12 and for varying values of r. We see that on 
increasing r, ^^^■^''~^(^) approaches a universal form until truncation effects take over and the 
TCSA approximation breaks down. 

There are two sources of error in the TCSA estimate ^^^''')(^). Firstly, TCSA gives 
the function ^^^'''HO x ((2//7r) sin(7r^//))~^/^ in equation Q as a Fourier series truncated 
at the 2A^*'^ term, which leads to the usual errors associated with truncation of Fourier 
series. Secondly, the coefficients /n^'^^ appearing in the truncated Fourier series are only 
approximately calculated. 
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-0.072820 


-0.072397 


-0.072162 


-0.072053 
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-0.015439 


-0.013388 


-0.003932 


-0.012698 


-0.012618 
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-0.0004456 


-0.004562 


-0.009036 


-0.003828 


-0.003715 
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-4.406-10-5 


-0.0001517 


-0.001821 


-0.001572 


-0.001524 


5 




-2.009-10-5 


-6.377-10-5 


-0.002017 


-0.0007604 


6 




-3.59-10-6 


-9.89-10-6 


-3.15-10-5 


-0.000480 
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-1.74-10-5 
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Table 

The TCSA coefficients for r = 8 and 4 < iV < 12. 
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n 




r 
2 




4 




6 




8 




10 




12 





-0 


780435 


-0 


759482 


-0 


728325 


-0 


694979 


-0 


663779 


-0 


635945 


1 


-0 


0049452 


-0 


0226611 


-0 


0476133 


-0 


0720535 


-0 


0923617 


-0 


108094 


2 


-0 


0005305 


-0 


0027410 


-0 


0068911 


-0 


0126176 


-0 


/A ^ /A£AiAi^/^ 

0192286 


-0 


026052 
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0001405 
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0007379 
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0019225 
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0037146 
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0060376 
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0087501 
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5 


-2 


796-10-5 


-0 


0001469 


-0 


0003855 


-0 


0007604 


-0 


0012814 


-0 


001952 


6 


-1 


642-10-5 


-8 


691-10-5 


-0 


0002332 


-0 


0004801 


-0 


0008648 


-0 


001443 


7 


-3 


10-10-8 


-8 


01-10-^ 


-4 


98-10-^ 


-1 


73-10-5 


-4 


44-10-5 


-9 


48-10-5 


8 


-4 


13-10-9 


-1 


14-10-^ 


-7 


85-10-'^ 


-3 


06-10-6 


-8 


76-10-f^ 


-2 


07-10-5 


9 


-1 


15-10-9 


-3 


20-10-8 


-2 


24-10-^^ 


-8 


95-10-^ 


-2 


64-10-*^ 


-6 


45-10-6 


10 


-4 


01-10-1° 


-1 


12-10-8 


-7 


88-10-8 


-3 


20-10"^ 


-9 


66-10-^ 


-2 


44-10-6 


11 


-1 


44.10-10 


-4 


02-10-9 


-2 


85-10-8 


-1 


17-10-^ 


-3 


64-10-^ 


-9 


50-10-^ 


12 


-4 


30-10-" 


-1 


22-10-9 


-8 


93-10-9 


-3 


84-10-8 


-1 


26-10-^ 


-3 


54-10-^ 


sum 


-0.786151 


-0.786151 


-0.786151 


-0.786151 


-0.786151 


-0.786151 



Table |: 

The TCSA coefficients /^^'''^ for = 12 and 2 < r < 12. 
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Figure [T| 

Plots of G{0 and X(^) against ^ from TCSA truncated to level 12. The upper set of lines are 
the (symmetric) functions Q{C)^ and the lower set of lines are the (asymmetric) functions I(^) 
for h = —1/2, = 0. These are plotted for r — 2"/^ with —2 < n < 8. As r increases, the 
functions Q (^) and 1(f) both approach universal functions which we take to be the expectation 
values of M'^/^Lp{x) on a half-line with boundary conditions 1 and $(0) respectively. For large 
r, truncation effects start to intervene, as can be seen in the slight 'ripple' discernible for n — S. 
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3.2.2 The strip with {^{h),l) boundary conditions 



This case is only slightly more complicated. The unperturbed expectation value is again given 
in (|2.4)- The TCSA estimation of the expectation value in the massive model with boundary 
perturbations then takes the general form 

Af2/5(y{x)>,„^,„|^ (3.9) 

Tl — j — 1 

where h is related to the reflection factor parameter b by h = h{b), and where the functions 
9nj '^^ (/i) and /i^^'^^ {h) have to be evaluated numerically. Since the ground state | ) of the 
TCSA Hamiltonian lies in the h = —1/5 representation, the matrix elements of ip are given in 
terms of the two chiral blocks /i,/2; furthermore, since this representation has a null state 
at level 2, one can eliminate all states with more than one mode which leads to the fact 
that one does not need to use higher than the first derivative of the chiral blocks in ( p.9| ). 
Since there are rather many coefficients g^'^\ h^n'^\ we shall not give any explicit examples. 

In figure |l| we also plot T{^) against ^ for various values of / between 0.5 and 16, for 
the fixed value h = 0, 6 = — 1/2. The excellent agreement between and G{S,) that can 
be seen on the half of the strip r/2 < ^ < r for the larger values of r is a good sign that 
the TCSA estimates of the functions are converging to their correct values. Table ^ below 
includes results on the convergence in N of J(^)(^) for values of ^ between 0.01 and 1.0 . 



3.2.3 The strip with {^{hi),^{hr)) boundary conditions 

The calculation of the expectation value of ip on the strip with two perturbed boundary 
conditions {^{hi),^{hr)) is in principle the same, except that the functional form is rather 
more involved. The unperturbed expectation value is again given in ([27 



where the functions fi{0) are the strip chiral blocks ( |2.3D . However, the massive perturbation 
introduces terms proportional to the other two chiral blocks and their derivatives, so that the 
TCSA estimation of the expectation value in the massive model with boundary perturbations 
takes the general form 

H^^'^HO = m2/5(^(x))(^(,^),^(,^))|^ (3.10) 

B E E (iS^'c^. 'o a(?) + Co.. /;(f ) 

n=0 k=l 

where hi = h{bi), hr = h{br) and the functions jin'^\hi, hr) and k\^'^\hi,hr) are evaluated 
numerically. Again, since there are rather many coefficients, we shall not give any explicit 
examples. For large values of r, the two boundaries are essentially non-interacting - this 
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was already seen in figure |T|. We therefore see no new phenomena over those seen aheady; 
the TCSA estimates of ( M^/^ip{x) ) near the left boundary from the system with boundary 
conditions ^{hr)) are barely distinguishable from those from the system with boundary 

conditions (^(/i), 1). Further confirmation is contained in table ^, which presents data from 
the two situations. 

However, for small r, the presence of two perturbed boundaries can destabilise the vacuum 
even for values of the parameters hi and hr which are less negative than — |/icrit|) the value 
for which a single boundary destabilises the bulk vacuum on a half line For bi + br = 0, 
1 6; I < 2, the ground state and first excited state have an exact crossing at a finite value of r, 
while for bi + br > there is a finite range of r for which they become complex. We discuss 
this in section 3.2, where we give examples of the spectra, and plots of the boundary field 
expectation values {M^^^(j)), for systems with two perturbed boundaries. 





TCSA truncation level N 
4 8 




12 


16 


FF truncation level 
1 2 


n 

3 


0.01 


-1.19134 
-1.16788 


-1.16743 
-1.16261 


-1 

-1 


16362 
16234 


-1.16306 
-1.16239 


-1.61719 


-1.21289 


-1.14851 


0.03 


-1.35236 
-1.32317 


-1.32259 
-1.31662 


-1 

-1 


31786 
31627 


-1.31716 
-1.31633 


-1.60971 


-1.34554 


-1.31311 


0.1 


-1.48965 
-1.45280 


-1.45199 
-1.44456 


-1 
-1 


44604 
44414 


-1.44517 
-1.44421 


-1.58467 


-1.45438 


-1.44451 


0.3 


-1.52882 
-1.48563 


-1.48386 
-1.47625 


-1 
-1 


47722 
47586 


-1.47633 
-1.47592 


-1.52209 


-1.47810 


-1.47657 


1.0 


-1.40584 
-1.37786 


-1.37064 
-1.37380 


-1 
-1 


37158 
37424 


-1.3738 
-1.37447 


-1.37977 


-1.37549 


-1.37546 


Table 



















The TCSA and FF estimates of {(p{£,) )<s>{o), for varying distance ^ from the boundary from the 
Form-Factor approach truncated at particle number n and from TCSA truncated to level N with 
r = 12 on strips with b.c.'s (<i>(0), 1) (upper line) and ($(0),$(0)) (lower line). 
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3.3 The comparison of the TCSA and Form-Factor resuhs 

The Form-Factor method gives the expectation value of ^{x) on the half-plane. While we can 
think of the half-plane the infinite-width limit of a finite strip, this is not accessible directly 
using the TCSA method, which is limited to strips of finite width r and is expected to perform 
best near to r = 0. To enable a comparison of the two methods, we shall simply take the 
TCSA results on a strip of width r = 12, rather than extrapolating finite width TCSA results 
to infinite width. (The error for small values of ^ from taking TCSA results at r = 12 should be 
much less than the one-particle FF contribution in the middle of the strip, which is ~ 0.2%). 

In figure ^ we show results for two cases: the boundary conditions (1) and (<I>(0)). We 
give the Form- Factor expansion (with g = g/2) truncated to one-, two- and three- particle 
states, and the TCSA data from truncation to level 16 (with r = 12), and they are clearly in 
excellent agreement. We also show the Form- Factor expansion assuming g = g, and it is clear 
that this is wrong. In table ^ we give some numbers illustrating the convergence in TCSA 
truncation level and in form-factor truncation level for various values of .01 < ^ < 1. 




Figure H 

Comparisons of {M'^^^ip{x) ) on a half-plane with boundary condition (1) (upper lines) and 
($(0)) (lower lines). The points are the TCSA data, the dashed lines the FF result up to 1 
particle with g = g/2, the dot-dashed line the FF result up to 2 particles, the solid line the 
FF result up to 3 particles. The dotted lines are the FF results up to 3 particles with g = g. 

The boundary condition '^(0) (6 = —1/2) was chosen because it is for this value of b 
that the accuracy of the TCSA is highest; comparisons between FF and TCSA results for 
the further values b = —1 and 6 = can be found in figure 10, in section 6T below. The 
curve there for 6 = is particularly interesting, as this is the case, mentioned at the end of 
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section 2, for which the only difference between the and the 1 boundary states is the 

sign of the odd-particle-number contributions. 

The two- and three- particle form factor expressions in figure |2| are barely distinguishable 
from the TCSA data, and so in figures ^ and ^ we plot log(( M'^I^Lp{x) )) against log(^) for 
the 1 and $(0) boundary conditions. We also show the leading behaviour 

(m2/V(x))$(0) = (2e)i/^<X('^)*{0) + (2O'/''X(l)'f(0) + 0(^''/') , (3.11) 
(m2/V(x))i = (20'/'<X(t)i + 0(C^'/'), (3.12) 

of small-^ expansions obtained from a perturbative treatment of the structure functions. 
(We intend to report on this approach elsewhere |Q.) We see that the three-particle FF 
approximation already agrees very well with these expansions for — 3 < log.^ ^ ~2. 
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Figure |3|: 

Plots of log|(M2/5(p(x))i| VS. log(^). The 
points are TCSA data, the dashed, 
dot-dashed and solid lines are from the FF 
expansion truncated to one, two and three 
particles resp., and the dotted line is (3.12) 
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Figure ||: 

Plots of log I ( M2/V(x) )*(o) I vs. log(0 
points are TCSA data, the dashed, 
dot-dashed and solid lines are from the FF 
expansion truncated to one, two and t hree 
particles resp., and the dotted line is (3.11 
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4 Excited states and energy density 



The TCSA and FF methods are not restricted to the ground state expectation values of the 
field ip{x,y). In this section we give a couple of examples of their wider applicability. 

4.1 The expectation value of ip in the first excited state 

It is just as easy in the TCSA to find the first excited eigenstate | 1 ) as to find the ground state 
and to find the corresponding expectation value of (p. The result is somewhat less accurate 
than in the ground state, and this accuracy decreases as higher and higher excited levels are 
considered. In figure ^ we show (as points) the TCSA result near the <&(0) boundary of the 
strip with boundary conditions (<5(0), 1) with r = 14 and truncation level 14. This state 
corresponds to the boundary bound state of energy ei = Mcos((6+ l)vr/6) = \/3M/2, and so 
we expect that near the $(0) boundary the expectation value is approximately given by the 
expectation value in the first excited state for the semi-infinite geometry. It turns out that 
this expectation value can also be obtained using FF techniques, with results that are shown 
in the various curves on the figure. The calculation relies on an idea of analytic continuation 
between states; similar methods were used to find TBA equations for excited states in Q. 



. 5 




-1.5 - 



Figure ||: 

Comparisons of {M^^^ip(x) )$(o) in the first excited state. The points are the TCSA data 
from truncation level 14 and r = 14. The dotted line is th e FF expression (3^) with 
b = 9/2; the dashed line is the correction by the terms in (4.2) from the poles at 
= ±(2— &)7r/6, and the soli d lin es are t he r esults of incl udin g all corrections by the 
addition of the terms up to (O), up to (|4.3| ) and up to (15) respectively. 



The fact that the refiection factors Rb{0) obey a curious continuation property was already 
observed in |8|. While the physical parameter h(b) is invariant under under b — > 4—6, the 
refiection factor is not; instead, it is interchanged with the refiection factor of the first 
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boundary bound state. This suggests that the expectation value of the field ip in the first 
excited state can be obtained by continuing the FF expression ( p.6| ) from the domain —3 < 
6 < 2 to the region 2 < 6 < 5. 

In figure ^ the result of the substitution of 6 = 9/2 in ( ^^ ) is shown as a dotted line. 
Clearly, it is a long way from the corresponding TCSA results for {lp)^^\ The explanation is 
simple: there are poles in Ki,{9) whose positions depend on 6, in particular at = ibi(2— 6)7r/6 
and 6 = ibi(4— 6)7r/6. These two pairs of poles cross the integration contour (the real axis) at 
b = 2 and 6 = 4 respectively, and contributions from both pairs must be added in explicitly 
to recover the correct analytic continuation in b of {(p) to b = 9/2, corresponding to h = 0. 

We can regard the contributions from these poles as directly affecting the exponential in 
( |3.4] ). We denote the positions of the 'active' poles (those which have crossed the integration 
contour during the continuation) by Oi(b), and the contribution to the contour integral from 
the pole in Ka[9) at 6i{b) by ki{b) - this will be ±1 times the relevant residue, depending on 
whether the contour was crossed from above or below when the pole became active. Then we 
can associate the following state in the full-line Hilbert space with the first excited state on 
the half-line: 



exp 



5.^(0) + - Y,h{b)A(-emm{b)) + j —K^{e)A[-e)A[e) 



o; 



(4.1) 

Expanding this out and inserting the appropriate form factors, the first corrections to the 
expectation value coming from the residue terms are 

^ J] A;,(6) (0195(0)1 -0i(6),e,(6))e-2M.co.h^^.(6) ^ ^4 2) 

i 

'-^Y^Hb) (0 1(^(0)1 0,-0,(6), 0,(6) )e-^^^(^+2coshe.(6)) ^ ^4 3) 



-lY.hm.ib) (0 1^(0)1 -0,(6), 0,(6), -0,-(6),%(6))e-2*^^(™^^^^'(^)+^°^^^^('')) , (4.4) 

oE^i(^) / ^^aW(0|v^(0)| -0,(6), 0,(6), -0,0)e-2*^-(-^i^^»('')+™^^^) , (4.5) 
^ ^ J-oo 47r 

where the last two terms can be seen as coming from poles in the four-particle contribution 
to the boundary state which hitherto we have neglected. Also, it turns out that the term 
(|4.4D always gives zero, due to the particular relative positions of the poles at 0i and 02. 
We have shown the result of correcting (3.4) by the dominant correction (the term in ([4. 2]) 



coming from the poles at = ±(2 — 6)7r/6) as a dashed line on figure |5|, and the result of 



adding all terms up to (O), up to ( [4. 3D and up to (4^) as solid lines. It is clear that these 



are converging rapidly to the TCSA value. 
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4.2 The expectation value of the energy density e 



As a second example, we use the TCSA to find the expectation values of the energy density 
e(x) on the strip for which 



H 



R 



e{x) dx , e(x) 



1 

2^ 



[T{x) + T{x)\ + \^{x) . 



One has to be rather careful about the specification of the operators in this expression. Here, 
we mean by T{x) and T{x) the 'bare TCSA' quantities - in other words, their expectation 
values are computed in any given state using the matrix elements of the CFT operators T 
and T between the eigenstates |n) of the perturbed Hamiltonian, themselves expanded in the 
basis of CFT states. This is the same procedure as was used for computations of ( ^{x) ) 
earlier, but some new issues arise in this case, to which we hope to return in [^. 

Leaving these questions to one side, the operator e allows us to see directly that 'boundary 
bound states' are indeed localised at the boundary. In figure |6| we plot the difference in the 
energy density between the first few excited states and the ground state for the system with 
boundary conditions (^>(0),1) for r = 8, calculated using TCSA truncated to levels 12, 14 
and 16. The result of truncating the Fourier-like expansions is very evident here, the TCSA 
estimates having distinct high-frequency ripples (varying with truncation level) superposed 
on the overall function. 

From the analysis in ||5|,^, for large r the first excited state is a boundary bound state, 
and the next several excited states are single-particle scattering states. This is borne out 
by figure ^, where we see very clearly that the first excited state corresponds to a particle 
trapped on the left ($(0)) boundary and decaying exponentially across the strip, while the 
higher excited modes are well spread across the strip, attracted to the $(0) boundary and 
repelled by the 1 boundary. 




Figure ||: Plots of M-2(( e(^) )„ - (e(C))o) vs. ^ for the first 
few excited states on the system with boundary conditions 
($(0), 1) and r — %. The solid line is the 1st excited state, 
the dash-dotted line the 2nd, and the dotted Une the 3rd. 
For each Une, TCSA results truncated to levels 12, 14 and 
16 are superimposed to give an idea of the errors. 




Figure |^: An attractive and 
a repulsive scattering process 
with the same classical time 
delay. 



It is impossible to decide whether a boundary {Ba) is repulsive or attractive purely given the 
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reflection factors Ra{0). The <J>(/i(0)) and 1 boundaries have the same reflection factor but 
from flgure ^, the $(/i(0)) boundary is attractive and the 1 boundary repulsive. This is the 
quantum analogue of the inability in a classical theory to determine whether a boundary is 
attractive or repulsive given only the time-delays - the two processes illustrated in figure ^ 
have the same time delay but clearly one describes attraction to the boundary and the other 
repulsion. 



5 The one-point functions of the boundary field 

We now turn to the one-point functions involving the boundary perturbing field (p. In prepa- 
ration for the main calculations, consider first the partition function Z^p of the model on 
cylinder of length R and circumference L, with boundary conditions a and (3 imposed at the 
two ends. Formally, we can write 

-Ably 



Zap = J[V^e-^^^^, (5.1) 

where J"[2?^'] implies that a functional integral over all bulk and boundary degrees of freedom 
is taken, and Ably denotes the combined bulk and boundary action. For the pair of boundary 
conditions {a,f3) = {^{hi),^{hr)), this can be written as 

pR pL pL pL 

Ably = Abcft + >^ dx dyip{x,y)+hi dy (f){0,y) + hr dy(f){R,y) (5.2) 
Jo Jo Jo Jo 

where Abcft is an action for the A^2,5 conformal field theory on the cylinder with conformal 
boundary condition <I> at the two ends. (For the other two pairs of boundary conditions the 
expression is similar, but lacks one or both perturbing boundary fields.) 

The behaviour of Z^p as a function of R and L is complicated, but if both are much larger 
than all bulk and boundary scales, then, up to exponentially-small corrections, 

log Zap ~ -RLS^uik - Lfa - Lffs (5.3) 

where fa and fp are the extensive parts of the boundary free energies, and i?buik the extensive 
part of the bulk free energy. For the scaling Lee- Yang model, Shuik = —M'^/{A\/^) [p4| . 

Given Z^p , the (normalised) one point functions of the field (p can be simply obtained by 
differentiation: 

(<Av.)cyl = -^^log^«,. (5.4) 

However, while it was shown in H] that the partition function was numerically accessible 
via the TCSA, at the current state of technology Zap (and hence {4')cyi) is not directly 
computable by means of the TEA. Contact with this 'exact' method can instead be made in 
certain limits, and these are best discussed using a Hamiltonian formulation. 

In fact, there are two alternative Hamiltonian descriptions of the partition function. In 
the so-called L-channel representation the role of time is taken by L: 

Zap = Tr^,^_^,e-^^<^''(^'«) (5.5) 
= ^M-LEf'^M,R)), (5.6) 

E„Gspec{H„0) 
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while in the R-channel representation the role of time is taken by R: 



Z. 



(a I exp(-i?i7circ(M,L))|/3 



espec{Hcirc) 



eM-RET{M,L)) . 



(5.7) 



In ( |5.7| ) we have used the boundary states \a) and |/?), and the eigenbasis {iV'n)} of f^drc , the 
Hamiltonian propagating states living on a circle of circumference L. The two decompositions 
are illustrated in figures ^ and 



L 

(time) 



R — 
(space) 



Figure ^: 

The L-channel decomposition: states \xn) live 
on the dotted line segment across the cylinder. 



L I 
(space) 



(time) 

Figure |[ 

The R-channel decomposition: states \ipn) live 
on the dotted circle around the cylinder. 



For the rest of this section we will focus on results for the boundary field (j)i on the left- 
hand end of the cylinder, but with a trivial relabelling it is clear that analogous results for 
can be found. 



5.1 L-channel decomposition 

First we consider the L-channel representation, depicted in figure ^. This will enable us to 
find exact formulae for the expectation value of (f) on strips of finite width. Introducing the 



eigenbasis {|Xn)} of ^{afi) we diff'erentiate inside the trace (|5.5D to find 



)cyl 



{Xn\(t>l\Xn 
^ (XnlXn) 



■exp(-L£;^*"P(M, R)). 



Comparing with ( |5.6| ) 



(Xn|</'/|Xr 



±.Er'^M,R), 



(5.8) 



(5.9) 



iXnlXn) 

where {XnlfpilXn) / {Xn\Xn) is the expectation value of the field (pi on the left boundary, taken 
in the n**^ excited state. 

It was shown explicitly in [S] that, at least for small n, En^^^ can be computed using 



generalisations of the boundary TBA equations of |16|. For the scaling Lee- Yang model this 
involves a non-linear integral equation for a single function £{9): 



e{e) = 2r cosh 9 -log \af3 (9) + 



S{9-9_p) 
3(9 -9p) 



IC*L{9) , 



(5.10) 
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and an associated set of equations for the (possibly empty) set {0p,9p} of so-called 'active' 
singularities (cf. [^,01): 

Here, r = MR as in earlier sections, L{e) = log(l+e-^(^)), f*g{e) = ^ f^^dO' f{9-e')g{9'), 
and, for the ($(/i(6;)), boundary conditions, 

KfsiO) = K,^ {e)K,^ i-9) , IC{9) = -z^ log 8(9) , (5.12) 



with Kh{9) defined in p.2| ), and S{9) the bulk S-matrix (p.7|). The number of active singular- 
ities depends on the particular energy level under consideration; for some pairs of boundary 
conditions on the strip it is nonzero even for the ground state |^], in contrast to the situation 
for the more familiar TBA equations for periodic boundary conditions. 

The solution to ( 5.10| ) for a given value of r determines the function Cn{r) 



6 12r — 
cjr) = ^ d9rcosh9L(9) + i ^(sinhep - sinh6l„) , (5.13) 

^ J-oo ^ „ 



in terms of which £'^*"''(M, R) is 



^P(M, R) = S^^i^R + 4 + ft,. - ^cn{r) , (5.14) 



and, rewriting (5.9) in terms of 



{Xn\4>l\Xn) _ ( d/t\ t,ip^ ^pstrip i 

{Xn\Xn) \dh) dbl " 7T\h,,,,\cOs{7T{bi + l/2)/5) dbl " ' ^ 



In the expression ( 5.14 ), fi,^ and /{,,, are i?-independent contributions to En^^'^iM, R) from 



the two boundaries, and <?buik is the bulk energy per unit length. If the equations for the 
particular state under consideration contain active singularities {9p, 9p} whose positions do not 
tend to zero as r — > oo, then there will be further i?- independent contributions to En^'^i^M, R) 
coming from the second term on the RHS of ( ^.13 ) - these will be described shortly. However 



for ground states such contributions are always absent, and so (as anticipated by the notation) 
fb^ and can be identified with the extensive parts of the boundary free energies as defined 
in (^1^). The exact values of these quantities were extracted from the (ground state) TBA 
equations in [Q, both for the ^{h{b)) and the 1 boundaries, and are 

fb = + sin f ) M , /i = fk=o . (5.16) 

States \xn) lying above the ground state |xo) will generally be separated from |xo) by a 
finite energy gap, even at large R. At the level of the TBA, this gap is seen in the presence 
of the active singularities, mentioned in the last paragraph, whose positions do not tend to 
zero as r ^ oo. These give Cn{r) a linear growth in r = MR, which via (|1|) yields an extra 



constant term in the large- i? asymptotic of En'^^^iM^R). Physically, the gap arises from two 
sources: the state \xn) niay contain a number, k{n) say, of bulk particles bouncing between 
the two edges of the strip, and in addition there may be a boundary bound state sitting at one 
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or both of the boundaries. Suppose that the possible boundary bound states for the left-hand 
boundary are indexed by = 0, 1, . . . A:max(&«) with k = the boundary ground state, and 
likewise for the right boundary, and that for the state \xn) the left and right boundaries are 
in states ki{n), kr{n) respectively. Then the general behaviour for i? ^ oo is as follows: 

'P(Af, R) ~ fbuiki? + Mk{n) + 4 + + ek,in){bi) + efc^(„)(6,) , (5.17) 

where 6^(6) is the energy of the k^^ boundary bound state of the <I>(/i(6)) boundary condition 
(with eo{b) = 0, since A; = is just the boundary ground state). 

Taking the limit i? ^ cxd of the ground state expectation value allows the boundary 
expectation values on the semi-infinite plane to be recovered. For the state xo ; ^(0) = 
ki{0) = kr{0) = 0, and substituting ( ^Tzl) into (|5ll ) using (|5.16D , we obtain the exact value 



of the dimensionless expectation value {M^^^cp) in a semi-infinite geometry: 

6|/icrit| cos(vr(6+l/2)/5) ^ ^ 

For now we take b to be restricted to the 'fundamental region' — 3 < 6 < 2, in which case the 
above formula is indeed correct for the ground state. Its interpretation as b moves outside 



this region will be described in section 6.1 below. 

Returning to finite values of R, the asymptotic ( |5.17| ) no longer suffices to obtain the 
expectation value of the boundary field, and one has rather to differentiate the exact formula 
( |5.14| ), so that ( |5.18| ) becomes 



^ ^'^^^ 6|/ie.it| cos(.(6, + l/2)/5) ' ^'-^'^ 

where ^c(r) = ^co(r) can be obtained by differentiating the full TBA equation ( ^.101 ) (a 
similar idea was applied to the case of periodic boundary conditions in [^^). Restricting, for 
simplicity, to the ground state energy in the region where no active singularities are present, 
the terms involving the 9p in ( |5.10| ) are absent, and 



where e{9) solves ( 5.10| ), while r]{9) = -^£{d) can be obtained from the linear integral equation 



obi l-l-exp(e) 

In this way the final estimates for r]{9) and £{9) have roughly the same accuracy. 

We also estimated ( 0) using the TCSA. For this, we used a strip geometry with specific 
boundary conditions (a, (3) on the two edges. We then calculated the dimensionless expec- 
tation value ( M^^^<f)i )a,p{f) , as a function of the strip width r for finite truncation level 
N. While {M'^/^(j)i)a^,3 can depend strongly on /? for small r, as r increases, this dependence 
decreases, and as r ^ oo it approaches the half-plane value ( M^/^ In table ^ we compare 
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the TCSA estimates for ( M-'^/^c/) )q,^^ for a = $(0) , /? = 1 and r small with the numerical so- 
lution of the ground state TBA equations ( |5.9| , ^. 10 , ^.11 ) and ( 5.20 , 5.21 ). A similar agreement 
was found for other pairs of boundary conditions. 
In table |5| we also give 

(MV50) for the same boundary conditions for larger values of r 
to show the convergence to the IR value, and include several plots of (M^/^c/)) for various 
boundary conditions in figure 14. 



log(r) 





-8 


-6 


-4 


-2 





TCSA 
TBA 


-0.19684202180 
-0.19684202181 


-0.29365378867 
-0.29365378869 


-0.4380797960 
-0.4380797961 


-0.65350591 
-0.65350592 


-0.969216 
-0.969217 



Table |: 

The TCSA estimates, with truncation level N = 14, of {M^/^(p)a^f3 for 
a = $(0) , /? = 1 in the small r region compared with the numerical 
prediction from the TBA 










1 


log(r) 
3/2 




2 




5/2 


6 


-0.969212 
-1.321108 


-1 
-1 


129791 
202314 


-1.166399 
-1.179004 


-1 
-1 


175226 
174604 


-1 
-1 


182737 
175250 


8 


-0.969214 
-1.321116 


-1 
-1 


129788 
202350 


-1.166291 
-1.179061 


-1 
-1 


174404 
174612 


-1 
-1 


178406 
174563 


N 10 


-0.969215 
-1.321120 


-1 
-1 


129795 
202370 


-1.166281 
-1.179099 


-1 
-1 


174229 
174643 


-1 
-1 


176402 
174357 


12 


-0.969216 
-1.321122 


-1 
-1 


129801 
202382 


-1.166281 
-1.179124 


-1 
-1 


174121 
174679 


-1 
-1 


175346 
174337 


14 


-0.969216 
-1.321124 


-1 
-1 


129806 
175194 


-1.166288 
-1.179142 


-1 
-1 


174104 
174709 


-1 
-1 


175196 
174342 



Table § 

The TCSA estimates of ( Af^/^0)<j,(o),i('^) (upper data) and 
( M-^/^i?!) )$(o),4>(o) ('^) (lower data). These can be compared with 
the 'exact' IR value of -1.17459499975... 
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5.2 R-channel decomposition 

We now turn to the R-channel representation, depicted in figure |9|. This will ultimately lead 
to expressions for expectation values when the boundary field is placed at one end of a semi- 
infinite cylinder. We shall use the notation of |^ and denote the boundary states | a ) and 
I /3) as I ^{hi) ) and | ) respectively. 



If R is taken to infinity in (5.7) with all other variables held fixed, then the contribution 



of the ground state l^/^o) = 1^) will dominate the spectral sum. Thus 

Za^ ~ A^f3{M, L) eM-RE^""{M, L)) , (5.22) 
where Eq^^{M, L) is the ground state energy of //circ and 

^.MM,L) = ^^^^^. (5.23) 

If we now let L grow as well and compare with (^.3|), we see that the inner products appearing 
in ( p.23| ) will, in general, contain a term corresponding to a boundary free energy per unit 
length: 

(TTTj^ ^ = + log(5.(M, L)) . (5.24) 

On the other hand, this linear term can be extracted from the small-L behaviour of the 
functions log((7a(M, L)), for which 'L-channel' TBA equations^ were proposed in |1(:]. This 



is explained in Q, where a precise match with the earlier result (|5.16| ) was found. The 
consistency between these two determinations of fi, is in some respects a mystery, since, from 
other results reported in |P , there are good reasons to doubt the ability of the L-channel TBA 
equations of ||l^ to describe the full variation of \og{ga{M, L)) as a function of ML. 
Returning to the R-channel decomposition of the full partition function, we have 

( >cyi = -T^ E ^(^OlV'n ) ^^"'^^/^\^^ xp(-i?Er(M, L)) (5.25) 
The following identification was made in the massless |j^] and massive cases: 

( $(/l(6))|V'n ) = Yn{i7T'-^){ IL|^„ ) , (5.26) 

where Yn{9) = exp(e„(0)), and en{0) is the solution of the excited-state TBA equation (with 
periodic boundary conditions) for the state \^n) (see fH,^). Taking the limit R ^ oo one 
deduces that 

which gives the 1-point expectation value of (p acting on the end of an infinite cylinder of 
circumference L in terms of the function Y. One check on this formula is easily made: from 
the large- -L limit of the TBA equation for Y, we have 

logy(ivr^)~-AfLsinf = -L(/fe-/i) {L ^ oo) (5.28) 



^the terminology is unfortunately, but unavoidably, a little confusing - these equations are called 'L-channel' 
because their derivation proceeds via the L-channel representation of the partition function. 
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using ( 5.16| ) in the second equality; differentiating, the expectation value on the upper half- 
plane quoted in (|5.18| ) is recovered. 

There is evidence (see 0) that the full set of "excited" Y^s can be obtained from Y via a 
process of analytic continuation in the bulk perturbing parameter A. It then seems reasonable 
to suppose that ( 5.27| ) can also be continued, leading to the following general relation: 

^"^'^ - imi^n) =~l[Tb) ^,logY.{^n—) . (5.29) 

(Alternatively, one can obtain this simply by differentiating ( 5.261 ).) 

It turns out that there is a further consequence of the R-channel decomposition. Note 
first that (K2^ can be used to write the partition function as 

= E ^ ^^ti^ll^nT ^ ^-(^^^)^»(^^^) ^M-RE^iM, L)) . (5.30) 



Now, the Y's satisfy the functional relation [27| 

Yr^{e + i^)Yn{9-i^) = l+Yn{e), (5.31) 

and this suggests the following identity 

Z^{h{b+2)),'S>{h{b-2)) = ■^1,1 + ^i,<i>(/i.(fe)) • (5.32) 



This should hold for all M, R and L; from the L-channel decomposition (5^), it is equivalent 
to the following relation between the spectra of models on strips of equal widths but different 
boundary conditions: 

{i^f 'P(Af, R)Uihib+2))Mhib~2)) = {ET'^iM, i?)}i,i U *P(M, R)U,^^hib)) ■ (5.33) 

Preliminary numerical work confirms this rather surprising identity, but as yet we do not 
have a good physical understanding of its origin. However, in section |6.2| below it will be 
used to formulate an exact conjecture concerning the regions of the (6, b') plane for which the 
model with {^{h{b)),^{h{b'))) boundary conditions develops a boundary-induced vacuum 
instability. 
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6 Applications 



In this section we apply some of the results obtained above to elaborate a few further aspects 
of the boundary scaling Lee- Yang model. We start with the boundary flows of the semi- 
infinite system, and then turn to the way that the boundary-induced vacuum instability of 
the model is affected when the system is confined to a finite strip. 

We first recall the way the parametrisations of the TCSA, TBA and FF calculations are 
related. The physical parameter describing the boundary is, of course, /i, which is related 
to b by 

h = - \hcrit\ sin(7r(6 + l/2)/5) . (6.1) 

This physical parameter h is periodic in b with period 10, and we choose the fundamental 
domain to be —3 < b < 2. Thus for real b, h is restricted to the range — |/icrit| <h< |/icrit|- 
To reach real values of /i > |/icrit| we can continue in h to values h = —3 -|- ib with h real; to 
reach real values of /i < — |/icrit| we can continue in b to values b = 2 + ib with b real. 



As we pointed out in section 4T, if we formally continue in b outside the fundamental 
region we have to be careful, as quantities may not continue naively - in particular the 
continuation of the 'fundamental' reflection factor to 2 < 6 < 5 in fact describes the reflection 
properties of the boundary with the addition of (the lowest) boundary bound state. 

6.1 Boundary flows on the semi-infinite system 

It is quite straightforward to see the boundary flow, at the level of the one-point functions, 
using the TCSA. In ||5| the spectra for the strip with boundary conditions ($(/i),l) were 
calculated for several values of h, and they are consistent with the idea that this spectrum 
is real for all real h > — |/icrit|- Thus we can calculate the expectation values of fix) and cj) 
on the boundary ^{h) by looking at the large r limit of calculations on the strip with bcs 
(<&(/i),lL). The only restriction is that the TCSA errors increase sharply with \h\, so that 
TCSA results are restricted to a small range of b values centred on 5 = —1/2. 

The FF and 'exact' (TBA) results have more interesting properties. Recall that the FF 
calculations are formally functions of 

tan((fe 2)7r/12) 
tan((5 - 2)7r/12) 



9<s>ib) = oL /.o; 91 , KbiO) = Rtiq - 9) . (6.2) 



The continuation to large positive value of h through b = —3 + ib works well, as one can 
easily verify that none of the 6-dependent poles in Kb(6) cross the integration contours in the 



FF integrals ( p. 6]) , and so none of the subtleties described in section above arise. The 
formulae 

Jim g^{-3 + ib) = <7i , Jim i^(_3+,^)(e) = Ko{9) = K^{0) , 

b^oo b~*oo 

can therefore be substituted directly into the form factor expansion, establishing the result 

hm (M2/V(x))$(,) = (M2/V(x))u (6.3) 

as an exact identity. One can similarly discuss the continuation of the expectation value of 
the boundary field, finding, as expected, 

Jim (mV5^)^(,(_3^^^))=0 = (mV5^)i. 
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To illustrate the result (O), in figure we plot ( M'^/^ip{x) ) close to the ^{h) boundary 
for various values of h. The FF results smoothly interpolate between the 1 and ^{h) bound- 
aries, and give good agreement with the TCSA results for small value of \h\. The h ^ oo 
limit is not directly accessible in TCSA, but we have extrapolated our results in h and we 
can see that (modulo an amplification of the Fourier-type truncation errors) it shows every 
sign of converging on the expectation value in the II be. 

As was discussed in for h > h{—l), in the massive system on the half-plane the 
boundary has no boundary bound states, and flows to the II boundary condition as h ^ +00. 
For < h < h{—l) this boundary has one bound state, and for — |/icrit| = h{2) < h < h{l) 
there are two boundary bound states; at h = — |/icrit| the ground state and first excited states 
become degenerate and for h less than this critical value — l/icriti, the system does not have a 
real vacuum. The presence of the bound states can be understood as particles being trapped 



near the boundary, as we have seen in section 4.2 



-0.5 



-1 . 5 




0.5 



1 . 5 



2 . 5 



3 . 5 



Figure 

Plots of ( M^^^ip{x) )<s>{h{b)) with boundary condition The solid lines 

are the Form-Factor results truncated at three particles with b taking values 
(from second line from top to bottom) —3 + 3i, —2, —1, —1/2, 0, 1/2. The points 
are the TCSA results !{£,) with truncation level 14 and r = 8. We have also 
included a plot of {AI^/^(p{x))ji (top line), the extrapolation of the TCSA data 
to = 00 (long dashed line) and the bulk value (short dashed line) for 
comparison. 



The continuation of the FF results beyond 6 = 2 was already discussed in section ^T 



The boundary expectation values can also be continued, and here we find a repetition of the 
reflection factor results - the ground state and first excited state expectation values swap 
under 6 — > 4 — 6, and the second excited state is invariant. To recall, the boundary free energy 
and the excitation energies of the two lowest lying states (for 1 < 6 < 2) are 

= (^^^ + sin(7r6/6)) M , ei(6) = M cos(7r(6 + l)/6) , 62(6) = M cos(7r(6 - l)/6) . 
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cos(7r 


6/6) 


6|/lcrit| / 


sin(7r(6 


- 2)/5) 


5 ^ 


cos(7r(4 


- &)/6) 
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sin(7r(6 


- 2)/5) 


5x/3 ^ 


sin(7r(6 


- 2)/6) 


6 /icrit j 


sin(7r(6 


- 2)/5) 



This gives (cf. eqs. ( ^.15 ) and ( 5.171 )) the expectation value {M^/^cj))^"''^ in these three lowest 
lying states as 

(MV50)(1) ^ MV5-d-(/, + ei) = 

Note that (</')(°) = (0)^^^ at the threshold for the first excited state 6= - 1, and ((/))(°) = 
{4>)^^^ at the threshold for the second excited state 6 = 1. 

Unfortunately the TCSA does not give very good results for the expectation values in the 
excited states for large r, and while extrapolations in r and truncation level indicate that 
these results are indeed correct, there seems little point in showing any plots. 

6.2 RG flows on the flnite size strip 

In 1^ the spectrum of the model with (1,$) boundary conditions was described in some 
detail, but results for the $) system were not presented. The analysis of these results 
provides a nice application of the spectral identity found in section |5.2| above, and in this 
section we describe how this goes. We begin with some results from the TCSA. 

In figures pT|-[T3| the finite size spectrum is plotted for the system on a strip with boundary 
conditions (<I>(/i(6;)), ^{h{br))) with bi = br = b taking values —1/2,0, 1/2. Observe that the 
ground and first excited states cross for 6 = 0, and for 6 = 1/2 an interval appears in which 
the ground state has left the real spectrum. However, so long as h > — |/icrit| the large-r 
spectrum stays real. 

In fi gure |14| we plot the log of the expectation value of the boundary field [M ' 4')a,p 
against log(r), for the two pairs of boundary conditions (a,/3) = (<I>(/i(6)), 1) and (a,/3) = 
($(/i(6)), $(/i(6))) for 6 taking the values —1/2,0,1/2. For large r, the expectation values 
in the two boundary conditions converge to the same value, as we expect; for large r the 
influence of the right boundary on the left boundary decreases and the expectation value 
tends to the half-plane value, {4>i)a,f3 — >r^oo {4')a, independent of /3 (provided that the 
system is not destabilised by the boundary condition (3). Conversely, for small r, {(p)a,i3 



tend to the conformal limits (2.5) which are governed by the UV fixed point of the boundary 
flow. On figure ^ the conformal limits are shown as dotted straight lines - the expectation 
value for (a, /3) = (<I>(/i(6)), 1) all converge to the lower straight line in the UV, and those for 
{a, (3) = ($(/i(6)), <I>(/i(6))) all converge to the upper line. 

The most interesting behaviour is that shown by ( M^/^</>)$(/j(b)),3>(h(b)) for intermediate 
values of r. For 6 = we saw in figure 12 that the ground state and first excited states cross 
at one point; at the same point {(f)) diverges with a characteristic 'A' behaviour. For 6=1/2 
the ground state drops out of the real spectrum for a finite range of r; (0) diverges at the 
edge points of this range. 
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Figure ^ The spectrum of the model 
on a strip with boundary conditions 
bi = br — —1/2, plotted against r. 




Figure ^ The spectrum of the model 
on a strip with boundary conditions 
bi = br = 1/2, plotted against r. 
The dashed lines indicate the real part 
of a pair of complex conjugate 
eigenvalues. 



Figure |l^: The spectrum of the model 
on a strip with boundary conditions 
bi=br = 0, plotted against r. 




Figure |T|: Plots of log(-( Afi/s^)^,,^)) 
against log r for pairs of boundary conditions 

and for 
& = — 1/2 (long dashed lines), 6 = (solid 
lines) and 6=1/2 (short dashed lines). Also 
shown (dott ed li nes) are the conformal 
expressions ([2.51). See text for more details. 
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The appearance of level crossing can be understood by making use of the spectral identity 
( 5.33| ). Shifting 6, and using the symmetry of the model about 6 = 2, the identity can be 



written as 

'P(M, R)U^hi^b))Mhm = {Et'^'iM, U {Et'^'^iM, R)}^Mh(2~b)) ■ (6.4) 

Consider now the two ground state energies which appear on the RHS of this equation. 
Depending on their relative values, one or other will correspond to ground state of the model 
on the LHS. Two limits are easily analysed. As i? ^ 0, the fact that the (1,$) conformal 
Hilbert space contains the field ip while the (1, 1) does not shows that 

R^O : En'^'^lniTi) ~ +oo 



-E^o*"'^l(i,*(/i(2-6))) ~ -OO (6.5) 



^0 1(1,1) 
l(l,$(h(2-6))) 

In the opposite, large- i?, limit, the behaviours follow from ( |5.17| ): 

^^oo : -^r'^ki.i) ~ £bnlkR+iV3-l)M , 

-£^o*"^l(i,<i>(/^(2-6))) ~ -fbuik-R + (VS - 1 + sin ^2^^)M . (6.6) 



In the latter formula, without loss of generality, we took b to be positive. Comparing ( |6.5D 
with (|6.6| ) shows that the relative values of the two ground states on the RHS of (^^) swap 
over when going from small to large values of R whenever b is less than 2. Since these states 
can be identified with the two lowest-lying levels of the (^>, $) model on the LHS of ( |6.4D , we 
see that a level-crossing in this model is inevitable in all such cases. Thus the point (0, 0) in 



the {bi,br) plane, shown in figure 12, belongs to a whole line of points (6, —b), \b\ < 2 which 
also exhibit a level crossing. The fact that the relevant two states are taken from the spectra 
of distinct models on the RHS of ( |6.4D prohibits their mixing and ensures that the crossing 



will be exact. Once the line bi + br = is left, the identity (6.4) can no longer be invoked 
and the exact level crossing is lost, as can be seen in figures |Tl| and 13. Observe in the first 
of these the lowest two levels remain real, while in the second there is an intermediate range 
of R for which their energies become complex. Physically this can be explained as follows: in 
the first situation, bi + b^ < and the boundary fields are less strong than in the 'marginal' 
case of bi + b^ = 0, and hence have less chance to destabilise the model; and in the second 
the story is reversed, the boundary fields are stronger, and there is therefore a possibility 
of a vacuum instability for some finite values of R. Once this has been understood it is 
reasonable to conjecture that for all points {bi, br) in the fundamental domain —3 < bi,br < 2 
with bi + br > the model with ($(/i(6/)), $(/i(6^.))) boundary conditions exhibits a vacuum 
instability for some range of system sizes, while for the points below this line, the spectrum 
remains entirely real at all values of R. 



This picture is confirmed by TCSA plots analogous to figures 11-13 taken at various 



other values of bi and br, and leads to the phase diagram for the model shown in figures 15 



and 16. The change in coordinates from {bi, br) to {hi, h,-) in passing between the two figures 



transforms the line segment (6, —b), \b\ < 2 into the portion of the ellipse {hi+hr) / sin + 
(hi—hr)'^/ cos^ = on figure |l^ which touches the (shaded) region of instability. (Along 



the rest of the ellipse, the spectral identity (|6.4[) also holds, but does not imply a level 



crossing.) Note also that whenever either hi < /icrit or hr < /icritj the vacuum is already 
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unstable in infinite volume and so this region can immediately be shaded in, without the 
need to appeal to more subtle arguments. No markedly new features emerge in the spectra 
for 6/ 7^ br, so we will not show any further plots. However, we note that the value of R at 
which the level-crossing occurs for models on the line bi + br = diverges as bi (or br) tends 
to 2. This is as one would expect, since in this limit the value of one of the boundary fields 
is approaching — l/icriti, and the corresponding boundary bound state is becoming degenerate 
with the vacuum in infinite volume. 




Figure |]i 

The fundamental region in the plane, 
showing (shaded) the region within which, for 
at least one value of the strip width, the model 
exhibits a boundary-induced instability. 




Figure [16|: 

As figure [is], but plotted on the {hi,hr) plane. 
The inner dashed line, a symmetrically-placed 
square of size 2/icrit, indicates the extent of the 
region covered for real values of bi and br ■ 



7 Conclusions 

In this paper we have given a rather detailed analysis of the one-point functions of both bulk 
and boundary fields in a simple but nonetheless non-trivial integrable quantum field theory, 
the scaling Lee- Yang model. A number of different techniques have been explored, with 
results which have been shown to be in good accord. Previous work on this topic has tended 
to use just one method (most usually the form- factor expansion) and thus the consistency 
that we have found with between the different approaches is an important confirmation that 
the earlier studies have been well-founded. For this purpose the use of the scaling Lee- Yang 
model as a testing-ground is very natural, but we feel that it would be very worthwhile to 
extend this work to further theories, and work on such matters is currently in progress. 

There are a number of open questions that arise. The modification to Ghoshal and 
Zamolodchikov's boundary state that we were forced to make in order to reconcile FF and 
TCSA results remains a numerical observation, for which we have no terribly compelling 
physical argument. Furthermore, in section 4.2 we noted that the the boundary states could 
sometimes be analytically continued, yielding 'excited' boundary states associated with the 
boundary bound states of the semi-infinite line. A systematic understanding of off-critical 
boundary states, both in infinite and finite geometries, is still lacking and would presumably 
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shed some light on the questions raised by our observations. It also might help to understand 
at a more fundamental level how the modifications to the boundary TBA equations of |jl6|, 
found in ||5| by an indirect method of analytic continuation, arise. 

It would also be worth investigating boundary-to-boundary correlation functions such as 



^{il^l)il^r))-{<Pl){cPr). (7.1) 



These are accessible by a simple generalisation of the techniques explained in section 6.1 above, 
and might be amenable to comparison with, for example, the results of lattice simulations. 
We also remark that the one-point functions at the end of an infinite cylinder calculated in 
section 6.2 are essentially the analogues, for the boundary fields of a semi-infinite system, 
of the finite-temperature expectation values discussed in |1^, 1^]. In the model we discussed 



these were relatively easy to obtain as the field under consideration was also the boundary 
perturbing operator. For more general fields in more complicated models it would presumably 
be necessary to develop the discussion more along the lines of the work in |l^. This seems 
to us an important open problem which should certainly be investigated further. 

Finally, we remark that the partition function identity discovered in section 6.2 merits fur- 
ther study. In some senses it can be considered as a first off-critical extension of the identities 
between (sums of) conformal partition functions that can be observed on examining the lists 
of such objects provided in, for example, |§. A physical understanding of why such identities 
should exist is lacking, even in the conformal cases, and perhaps the broader perspective 
provided by the off-critical results will help towards this end. This alone should motivate the 
extension our work on off-critical boundary integrable models to further examples. 
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